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Abstract 

It is known that exact traveling wave solutions exist for families of 
(n + l)-states stochastic one-dimensional non-equilibrium lattice models 
with open boundaries provided that some constraints on the reaction rates 
are fulfilled. These solutions describe the diffusive motion of a product 
shock or a domain wall with the dynamics of a simple biased random 
walker. The steady state of these systems can be written in terms of 
linear superposition of such shocks or domain walls. These steady states 
can also be expressed in a matrix product form. We show that in this 
case the associated quadratic algebra of the system has always a two- 
dimensional representation with a generic structure. A couple of examples 
for n — 1 and n — 2 cases are presented. 

On the macroscopic level nonlinear hydrodynamic equations, such as the Burg- 
ers equation or the Fisher equation, might exhibit shocks in some cases [HE]. 
On the microscopic level some stochastic non-equilibrium one-dimensional lat- 
tice models with open boundaries might also develop shocks provided that some 
constrained on the reaction rates are fulfilled [21 [H [S] . 

In [Hj exact traveling wave solutions have been obtained for three families of 
two-states driven- diffusive models in one-dimension with nearest-neighbors in- 
teractions. These three models are the Partially Asymmetric Simple Exclusion 
Process (PASEP) [3 [51 [3] , the Branching-Coalescing Random Walk (BCRW) 
and the Asymmetric Kawasaki-Glauber Process (AKGP). By taking a Bernoulli 
shock distribution as initial distribution, the shock distribution in these models 
evolves in time into a linear combination of similar distributions with different 
shock positions. The time evolution of the shock position in these models is 
similar to the dynamics of a random walker moving on a one-dimensional lat- 
tice with reflecting boundaries ^ . This property of the traveling wave solution 
is not limited to the two-states systems with nearest-neighbor interactions. 
Shocks in three-states lattice models with nearest-neighbor interactions [TUl HH 
[T2I [T3l [M] and also two-states lattice models with next nearest-neighbors inter- 
actions [15] have also been studied recently. As long as the time evolution of 
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the product shock measure is equivalent to that of a random walker on a dis- 
crete lattice with homogeneous hopping rates in the bulk and special reflection 
rates at the boundary, the steady state of the model can be written as a linear 
superposition of these product shock measures. 

The steady states of these systems, on the other hand, can be obtained using 
the Matrix Product Formalism (MPF) first introduced in [5] and then developed 
in [inj (for a recent review see [Uj). According to the MPF the steady state 
weight of a given configuration of a two-states stochastic non-equilibrium one- 
dimensional lattice models with open boundaries can be written as expectation 
value of product of non-commuting operators associated with different types of 
particles and vacancies in that configuration. These operators should satisfy an 
associative algebra. By using a matrix representation of this algebra or using 
the commutation relations of the operators one can in principle calculate the 
unnormalized steady state weights. Whether or not this algebra has a finite- or 
infinite-dimensional representation in not known in general; however, it is known 
that some quadratic algebras have finite-dimensional representations on special 
manifolds in their parameters space. For instance it has been shown that the 
quadratic algebra of PASEP has finite-dimensional representations given that 
one is restricted to some special values of the reaction rates through some con- 
straints [IHl m] . Interestingly, it has been shown that the constraints for the 
existence of a two-dimensional representation are exactly those necessary for a 
single product shock measure as an initial configuration to have a simple random 
walk dynamics on the lattice [5D]. It has also been shown that if one takes a 
product measure of — 1 consecutive product shocks, the shock positions have 
simple random walk dynamics provided that the constraints for the existence 
of an Af-dimensional representation for the quadratic algebra of the PASEP 
are fulfilled. To see this one should adopt a slightly different definition of the 
shock measure. According to this new definition the shock is a product mea- 
sure with density 1 at the shock positions and intermediate densities between 
these sites [2D]. Hence the A^-dimensional representations of the stationary al- 
gebra describe the stationary linear combination of such shock measures. The 
same is true for both the BCRW and the AKGP [6]. In fact it is known that 
the quadratic algebras associated with the BCRW and that of the AKGP have 
two-dimensional representations under some conditions on the reaction rates. 
Surprisingly, it has been shown that a single product measure in these systems 
will have a random walk dynamics under exactly the same constraints pTj. In 
[llj a three-states model with open boundary has been studied and it has been 
shown that the traveling shock solutions exist under some constraints on the 
reaction rates. Later a generalization of this model was introduced and stud- 
ied in The authors in [10! have shown that the quadratic algebra of this 
generalized model can be mapped into the quadratic algebra of the PASEP pro- 
vided that the constraints under which the traveling shock solutions exist, are 
fulfilled. This means that the quadratic algebra of this generalized model has 
two dimensional representation under the same constraints. 
Our major motivation in this paper is to investigate how the steady states of 
one-dimensional driven-diffusive systems with open boundaries when written 
in terms of linear superposition of Bernoulli shocks are related to the matrix 
product steady states. In this direction, we will consider those {n + l)-states 
stochastic non-equilibrium one-dimensional lattice models with open bound- 
aries and nearest-neighbors interactions in which the time evolution equation of 
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a single product shock measure is similar to the evolution equation for a biased 
single-particle random walk moving on a finite lattice with reflecting bound- 
aries. Since the dynamics of the shock position, as a single-particle excitation, 
is simply a random walk it is easy to find the steady state probability distribu- 
tion function of these systems in terms of linear superposition of single shocks. 
We will then show that the matrix product steady state of these systems can 
always be obtained using two-dimensional representations of their quadratic al- 
gebras. These matrix representations have always a generic form regardless of 
the details of the microscopic reactions. In order to confirm our assertion we 
will give a couple of example for both n — 1 and n — 2 cases. 
In what follows we will first review very quickly the dynamics of a biased single- 
particle random walk on a discrete lattice of finite length and show how the 
steady states probability distribution of this system can be obtained. The 
time evolution of the probability distribution function of any configuration of a 
Markovian interacting particle system \P{t)) is governed by a master equation 
which can be written as a Schrodinger like equation in imaginary time 



dt 



\P{t)) = H\Pit)) 



(1) 



in which _ff is a stochastic Hamiltonian [8j. The matrix elements of the Hamil- 
tonian are the transition rates between different configurations. For the one- 
dimensional systems defined on a lattice of length L with nearest neighbors 
interactions the Hamiltonian H has the following general form 



L-l 



H = ^ hk,k+i + hi + Hl- 



(2) 



fe=i 



in which 



hk,k+l 

hi 
hL 



(3) 
(4) 
(5) 



For an (n + l)-states one-dimensional driven diffusive system, in which these 
states are associated with n different types of particles and empty sites, I is an 
(n -I- 1) X (n -|- 1) identity matrix and h is an (n -|- 1)^ x {n+ 1)^ matrix for the 
bulk interactions. The two matrices /i^''^ and ft.*^'-* are both (rt + 1) x (?i + 1) 
square matrices and determine the interactions at the boundaries. 
Let us assume that at f = the probability distribution function of the an 
(n -|- l)-states one-dimensional driven diffusive system is given by a product 
shock measure defined on a lattice of length L as follows 



|fc) = 



p^' 



\ 



rsL-k 



0< k< L 



(6) 



in which the density of particles of type i (i = 1, • • • , n) on the left (right) hand 
side of the shock position is p^*-* (pj'^ ) . The density of the empty sites on the 
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left and right hand sides of the shock position are given by p^^-* = 1 — J27=i Pi' 
and p"2^ = 1 — X]r=i respectively. As can be seen every species of particles 
has a shock profile and the centers of shocks lie on each other. We assume 
that under some constraints on the reaction rates of the system the dynamics 
of this state \k), which is given by ([T]), is similar to the dynamics of a biased 
single-particle random walker on a discrete lattice with reflecting boundaries. 
In this case, regardless of particle type, the shock position hops to the left and 
right with the rates 5i and (5,. respectively. The shock position also reflects from 
the boundaries. It reflects from the left boundary with the rate 5^ and from the 
right boundary with the rate 5i . The evolution of the state | k) is then given by 

H\k) = 6i\k-l)+5r\k + l)-{5i+5r)\k),{)<k<L (7) 
H\Q) = 6r\l)-6r\0), (8) 
H\L) = 6i\L-l)-Si\L). (9) 

Now one can construct a linear superposition of these single shocks to make the 
steady state of the system \P*) as follows 

L 

\P*)=Y,c,\k) (10) 

fe=0 

so that it satisfies 

H\P*)=0. (11) 
Using ([7)l- pT|) one can easily find the coefhcients c^'s as 

Ck = ^TTi ^-^^^ ,0<fc<L(12) 

^ ^ {Sr/6r) 

The equation ((TU]) gives the steady state probability distribution function of 
the system in terms of superposition of product shock measures. Because of the 
uniqueness of the steady state of our system, one should be able to find the same 
distribution function using other methods such as the MPF. In what follows we 
will briefly review the basic concepts of the MPF for an {n + l)-states system. 
According to the MPF the stationary probability distribution \P*) for a system 
of length L is assumed to be of the form 



( Do 



\V) (15) 



V Dn J 



in which Z is a normalization factor and can be obtained easily using the nor- 
malization condition to be 

n 

Z={W\{J2d,)^\V) (16) 

i=0 
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The operator Dq stands for the presence of a vacancy and the operator Di 
{i — 1, ■ ■ ■ ,n) stands for the presence of a particle of type i at each lattice site. 
These operators besides the vectors {W\ and \V) should satisfy the following 
algebra 



{W\ /i« 



( ^" \ 



\ Dn J 



V Dn J 



V Dn J 



V Dn j 



V Dn J 







( Do\ 








V Dn ) 




\ Dn J 



\v) = 



V Dn J 



\V) 



(17) 

in which Z?i's are auxiliary operators [111 HZ]- These relations define a quadratic 
algebra. Using the algebra or one of its matrix representations, one can calculate 
the probability of any configuration in the steady state. It is also possible to 
calculate the mean values of physical quantities such as the current of particles 
of different species or the correlation functions in the steady state. 
As we mentioned above, by requiring some constraints on the reaction rates 
of the system the dynamics of a single shock can simply be given by a single 
random walk dynamics. From there one can construct the steady state of the 
system. We should note that the results obtained from the random walk picture 
should be equal to those obtained from the MPF and in fact we should have 



:{W\ 



V Dn J 



\V)=Y.c,\k). 



(18) 



Let us assume that we have a matrix representation for the operators Di's and 
the vectors \V) and {W\. Finding the bra (fc'| which is orthonormal to \k) i.e. 
{k'\k) = dk'k helps us calculate the coefficients Ck- By defining 



\k') 



f ^0 \ 

Xi 
\ Xn J 



/ 2/0 \ 

2/1 

V Vn J 



(19) 



and requiring the orthonormality condition one finds that only four of 2n + 2 
parameters {xq, ■ ■ ■ , XmUo, ■ ■ ■ ,yn} are independent. One can simply find a 
solution by taking xq 



P2 



Xi 



2/1 



Vn 



in which pi — X]"=i Pi^ arid p2 = ^27=1 P2' the total 
density of particles on the left and the right hand sides of the shock position 
respectively. By multiplying {k'\ in (fT8|) from the left one finds 



-(1-P2) 
P2-P1 ' 



yo 



-pi 



and 



(0 



P2-P1 



Ck 



l{W\G'',Gt''\V) 



(20) 
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in which 

n n 

Gi=J2 ^» A and G2 = ^ y^D,. (21) 

1=0 i=0 

We have found that the foUowing two-dimensional representation 







do 




A=f ) . = l,---,n , (22) 



in which da = — X]"=i uniquely generates the same coefficients for k 
0, • • • , L as in HSll-dill) provided that 



do— = (P2 - Pi) — (23) 

l-(t)+(^) 

These two equations tell us how the reflection rates Sr and Si are related to the 
components of the two vectors \V) and {W\. One should note that in ([2^ . d^'s 
(i = 1, • • ■ , n) are free parameters. Since the matrix X)"=o diagonal in this 
representation, the normalization factor Z in can easily be calculated as 

Z = wiVi + (^)^w;2W2. (25) 
The matrix representations for the two matrices Gi and G2 are now given by 

do Sr_ \ 1^2 - \ -da 
P2-P1 Si / \ P2—P1 



Gi - ( do 5, ) , G2 - ( -do Q ) (26) 



with these properties 



- (|)'"'Gi , G^ - G2. (27) 



In summary whenever the quadratic algebra of an (n + l)-states system has 
a two-dimensional representation of the form (|22p one can calculate the coeffi- 
cients Cfe using PH]) and conclude that the steady state of the system is a linear 
superposition of single product shock measures with random walk dynamics. In 
order to confirm our assertion we will give a couple of examples. 
Let us first study the shocks in two-states systems which can be interpreted as 
particles (denoted by 1) and vacancies (denoted by 0) at each lattice site. The 
stochastic dynamics are defined in terms of transition rates. In this case the 
process is fully defined by the 12 rates Wij (i 7^ j) where i, j = 1, • ■ • , 4. In the 
basis (00, 01, 10, 11) the reactions in the bulk of the system are as follows 

01 ^ 10 W32, 1023 
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11 10, 01 W34, 'W2i 

10, 01 11 W43, W42 
10,01 -> 00 u;i3, W12 

00^10,01 ^31, ^^21 

11 ^ 00 U'14, W41. 

For injection and extraction of particles at the boundaries we introduce the four 
rates a, (3, 7 and 6. The reactions at the left and the right boundaries are 

0^1 ct, 7 and 
0^1 S, (3 

respectively. As the first example we consider the PASEP. It is well known that 
the steady state of the PASEP can be written in terms of superposition of shocks 
with random walk dynamics provided that some constraints on the reaction and 
the boundaries rates are fulfilled. In the bulk of the lattice the particles diffuse 
to the left and right with the rates W32 = x and W2^ — 1 respectively. For our 
continence we consider the injection and extraction rates for the left boundary 
as (1 — x)a and (1 — x)"f and also for right boundary as (1 — x)5 and (1 — x)(5. 
The associated quadratic algebra of this model can now be written as 

DiDo - xDoDi - e(l - x){Do + Di) 
(iy|(ai?o-7^i-0 =0 (28) 
{l3D^-5Do-S,)\V) = Q. 

It is assumed that Dq — ^{l—x)I and Di = — ^(1— a;)X in which! is a 2 x 2 iden- 
tity matrix. One can easily check the following two-dimensional representation 
satisfies ^ 



provided that 



/Ol(l - P2) &r _ /02(1 - P2) 



(30) 



P2(l-pi) ' 5i pi(l-pi) 

and f = ^2(1 — /32)- The boundary rates together with pi, p2 and x should also 
satisfy the following conditions 

"(1 - Pi) - 7Pi = Pi(l - Pi)(l - a;) 

(3p2-5{l- P2) = P2{1- P2){l-x). ^■''> 

In this case none of the elements ui, V2, wi and W2 in the vectors \V) and 
{W\ is zero. Finite-dimensional representations of the PASEP have widely been 
studied in [THl HH] . It is known that the PASEP has also an A^-dimensional 
representation provided that 

a;i-^ = K+(/3,<5)^+(a,7) (32) 

in which 

-u + V + I + ^ [u - V - ly + Auv 

K+{u,v) = . (33) 

2,u 



7 



Pi 

P2 



PJV 



N -1 



Figure 1: Sketch of a multiple shock structure in the PASEP with open bound- 
aries. The shock positions are at the sites 1, 2, • ■ •, and — 1. 



If one takes an initial state consisting of — 1 consecutive shocks (see Figure 
1) with the densities pi, p2, • • ■, and pjv, it turns out that the shock positions 



have simple random walk dynamics provided that the 



1, ■ • • , A^ - 1 and that 



pi+i(i-Pi) _ 
Pi(i-Pi+i) 



Pi 



Pn 



1 + K+(a,7)' 



^ for 

(34) 
(35) 



It is interesting to note that by choosing an appropriate form for the auxiliary 
operators Dq and Di an A^-dimensional representation of the PASEP 's quadratic 
algebra can be written as the following form 



(1 - pn) 



Do 



-(1 - PiV-l) 



(1 - P2) 



1 



■5,.,! 



in a basis in which Di is diagonal but not the Dq. In this matrix representation 
the quantities of dr,i and Si^i for i — 1, • • • , A^ — 1 arc the hopping rates of the 
ith shock position to the right and to the left respectively. This can easily be 
done by using a similarity transformation. It is worth mentioning that this 
representation is exactly the one introduced in [I9j however, it is written here 
in terms of the densities of the shocks pi, • ■ • , pjv and their hopping rates Sr^i 
and 6i^i for i = 1, ■ ■ ■ , N — 1. It can be easily verified that by taking N = 2 and 
applying a similarity transformation the two-dimensional representation (|22p 
can be recovered. One should note that the eigenvalues of the matrices are not 
changed under similarity transformations. 

As another example, we consider the BCRW as a two-states system with open 
boundaries and the following non-vanishing rates 



W34, W24, W42, W43, W32, W23, 



a, 



7, P 



(36) 
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and 6 equal to zero. The quadratic algebra of the BCRW is 



DoDo - DoDo = 

UJ23D1D0 + uj2iDl - {UJ32 + ^i2)DoDi = D^Di - DoDi_ 

-(W23 + UJ43)DiDo + + W32-D0-D1 ^ DiDo - DiDq 

uJisDiDo - (w24 + L^34:)Dl + LU42D0D1 = DiDi - DiDi 
{W\{aDo -jDi) ^ {W\Do ^ ~{W\Di 
PD,\V)=Do\V) = -Di\V). 



(37) 



It is known that for the BCRW the dynamics of a single product shock measure 
under the Hamiltonian of the system is a random walk provided that [B] 



-Pl C.J24+"34 



Pi 1^^42+1^43 

^ = ^ (38) 

7 = ^7^"" + (1 - Pl)uj32 - ^7f^W43 + P1W34. 

In this case the density of the particles at the right hand side of the shock 
position is zero p2 = 0. One can easily check that the following two-dimensional 
representation 

..^(;.;_,,)-^(-U;) (») 

with 5r = and 61 — (1 — pi)?i'32 + P1W34 satisfies ([57]) provided that the 



constraints (|38p are satisfied. 

For the AKGP, the non- vanishing parameters are W12, wia, ^42, W43, tL'32, a and 
p. The particles are allowed to enter the system only from the first site with the 
rate a and leave it from the last site of the lattice with the rate /?. In this case 
we have pi = 1 and p2 = 0. There are no additional constraints on the rates for 
this model. The dynamics of shock measure generated by the Hamiltonian of 
the system will be a simple random walk on the lattice and the shock position 
hopping rates are Si = a'13 and 5r = 0J43 [B]. The quadratic algebra of the 
AKGP is given by 

UJ13D1D0 + UJ12D0D1 = DqDq - DqDq 
-{lUi2 + UJ32 + ^^42)DoDi = DqDi - DqDi 

-(Wi3 + UJ43)DiDq + UJ32P0D1 = DiDq - DiDo , , 

U43D1D0 + ojMDi = DiDi - DiDi ^ ^> 

{W\aDo = {W\Do = -{W\Di 
PD,\V) = Do\V) = -Di\V) 

which has a two-dimensional representation given by 

One can also apply our procedure to the systems with n > 2. In what follows we 
investigate an open system with three-states at each lattice site associated with 
two different types of particles (denoted by 1 and 2) and vacancies (denoted by 
0). The particles of different types are allowed to enter and leave the lattice from 
the boundaries. There is also a probability for changing the particle type at the 
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boundaries. A couple of models of this type have already been studied in [inilTTl 
[13 [mill] and shown that a traveling wave solution might develop in the system. 
In this case the stochastic dynamics are defined in terms of 72 transition rates 
Wij {i^j) where i,j = 1, - • • ,9. In the basis (00,01,02,10, 11,12,20,21,22) we 
will consider the case studied in [TT in which the nonzero transition rates are 

W24, W42, ^37, W73 , wge , wgs , Wei, wie, Wis, wsi- (42) 

For injection and extraction of particles at the left boundary we introduce the 
rates : 

1^0 ai, 71, 

1^2 a2, 72, 

0^2 aa, 73, (43) 

and for the right boundary 

1^0 <5i, /3i, 

1^2 52, /?2, 

2 (53, /33. (44) 

The quadratic algebra of the system in this case is given by 

-(w6i + wsi)DoDo + wieDiD2 + wisDjDi = DqDq - DqDq 
-W42D0D1 + W24D1D0 = DpDi - DqDi 

-WnDoD2 + W37D2Do =_L'o^2 - Do_D2 

W42D0D1 - W24:DiDo = DiDo - DiDo 

wqiDqDq - (wi6 + WiiQ)DiD2 + wesDjDi = D1D2 - D1D2 

W73DQD2 - W3JD2D0 = D2DQ ~ D2D0 _ _ , . 

wjuDqDo + WS6D1P2 - {wis + wes)D2Di = D2D1 - D2D1 

D2D2 - D2D2 = DiDi - DiDi = _ 

{W\{-fiDo - {ai + a2)Di + 72 £'2 + D_i) ^ Q 

(I^|(a3i:'o + a2Di - (72 + 73)£'2 + i:'2) = 

(/3iDo - {5i + 52)Di + 132D2 - Di)\V) = 

{S^Da + 52D1 - {[32 + /33)^2 - D2)\V) = 

It has been shown that the shocks of the form © can develop in the system 
provided that some constraints are fulfilled. The shock in this model only hops 
to the left i.e. (5/^0 and &r = 0. In the following we show that the steady state 
of the system can be written as a matrix product state using two-dimensional 
representations of the quadratic algebra (|l5|) . We will study two different cases. 



In the first case we assume W42 — wis — wes — and also pj^'' — l,p'2^ — 0. 
The two densities p"^^ and p^^^ are free parameters. The matrix representation 
of the algebra (|l5l) is now given by 

in which do + '^1 +<^2 = 0. As can be seen the two operators Dq and D2 commute 
with each other but not necessarily with Di. The auxiliary operators Z?o, Di 
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and D2 are taken to be zero. The matrix element W2 in (j22p is zero in this case 
and the foUowing constraints have to be fulfilled 

5i = 82 = OLi ^ a2 ^ 0. (47) 

The rest of the boundary rates besides di and 6,2 should satisfy two linear 
equations 

/3idi + (/3i-/32)d2 =0 (48) 

(/3i + 53)di + (/3i + /33 + ^3)^2 = 0. (49) 

Assuming di,d2 ^ then those boundary rates which satisfy /3i/33 + /3i/32 + 
/32<53 = generate a solution. 

In the second case we consider p'^^ — pl,^^ — 0, p"2^ — 1 and p"^'' is a free 
parameter. The steady state of the system in this case is given by following 
two-dimensional matrices 

which is similar to (|46| except that d2 = Q- This means that the probability of 
finding second-class particles is zero in the steady state. Taking all the auxiliary 
operators -Dq, Di and D2 equal to zero, one can simply check that 150]) satisfy 
the algebra provided that 1^42 ~ and 

5i = 82 = ai ^ a2 ^ Pi ^ 5z = 0. (51) 

In this paper we have introduced a general procedure that can be used to inves- 
tigate those one-dimensional multi-species system with open boundaries whose 
steady states can be written as a superposition of single product shock mea- 
sures of the form ^ with random walk dynamics. Instead of applying the 
Hamiltonian of such systems on a single product shock measure © as an initial 
configuration and looking for the conditions under which this initial configura- 
tion evolves into a linear combination of product shock measures with different 
shock positions, as it is done in [TTl [131 El O'^s '^^'^ simply follow the 
following steps: 

1. Using the standard MPF we find the quadratic algebra associated with 
the one-dimensional open boundaries system. 

2. We investigate whether or not it has two-dimensional matrix representa- 
tions and that if these representations have the same structure as (P^ . 

3. If the answer to these questions is positive, we then calculate the coeffi- 
cients Cfc's in pn)) using ((^ . 

We should emphasis that one cannot use any arbitrary matrix representation to 
generate the coefficients Cfc's using (PUj) . Only when the two-dimensional matrix 
representation of the quadratic algebra is of the form ((^H) the coefficients c^'s 
have the right structure IT ^ - ITIl) . 

It would be interesting to investigate if this procedure also works for the systems 
with long range interactions |15| or even for the systems in which multiple 
product shock measures might evolve in them. As we have shown the generic 
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form of the matrix representation still holds for the PASEP with multiple 
shocks; however, it is not clear if it is true for other systems in general. It is 
also interesting to investigate if this property holds for the models defined a 
ring geometry or in the presence of second class particles. This is under our 
investigations. 

F. H. J. would like to thank V. Rittenberg, G. M. Schiitz and R. B. Stinchcombe 
for useful discussions and comments. 
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